4 results
Stokes flow singularity at the junction between impermeable and porous walls
- Ludwig C. Nitsche, Prashanth Parthasarathi
-
- Journal:
- Journal of Fluid Mechanics / Volume 713 / 25 December 2012
- Published online by Cambridge University Press:
- 24 October 2012, pp. 183-215
-
- Article
- Export citation
-
For two-dimensional, creeping flow in a half-plane, we consider the singularity that arises at an abrupt transition in permeability from zero to a finite value along the wall, where the pressure is coupled to the seepage flux by Darcy’s law. This problem represents the junction between the impermeable wall of the inflow section and the porous membrane further downstream in a spiral-wound desalination module. On a macroscopic, outer length scale the singularity appears like a jump discontinuity in normal velocity, characterized by a non-integrable $1/ r$ divergence of the pressure. This far-field solution is imposed as the boundary condition along a semicircular arc of dimensionless radius 30 (referred to the microscopic, inner length scale). A preliminary numerical solution (using a least-squares variant of the method of fundamental solutions) indicates a continuous normal velocity along the wall coupled with a weaker $1/ \sqrt{r} $ singularity in the pressure. However, inconsistencies in the numerically imposed outer boundary condition indicate a very slow radial decay. We undertake asymptotic analysis to: (i) understand the radial decay behaviour; and (ii) find a more accurate far-field solution to impose as the outer boundary condition. Similarity solutions (involving a stream function that varies like some power of $r$) are insufficient to satisfy all boundary conditions along the wall, so we generalize these by introducing linear and quadratic terms in $\log r$. By iterating on the wall boundary conditions (analogous to the method of reflections), the outer asymptotic series is developed through second order. We then use a hybrid computational scheme in which the numerics are iteratively patched to the outer asymptotics, thereby determining two free coefficients in the latter. We also derive an inner asymptotic series and fit its free coefficient to the numerics at $r= 0. 01$. This enables evaluation of the singular flow field in the limit as $r\ensuremath{\rightarrow} 0$. Finally, a uniformly valid fit is obtained with analytical formulas. The singular flow field for a solid–porous abutment and the general Stokes flow solutions obtained in the asymptotic analysis are programmed in Fortran for future use as local basis functions in computational schemes. Numerics are required for the intermediate-$r$ regime because the inner and outer asymptotic expansions do not extend far enough toward each other to enable rigorous asymptotic matching. The logarithmic correction terms explain why the leading far-field solution (used in the preliminary numerics) was insufficient even at very large distances.
Nonlinear drift interactions between fluctuating colloidal particles: oscillatory and stochastic motions
- E. J. Hinch, Ludwig C. Nitsche
-
- Journal:
- Journal of Fluid Mechanics / Volume 256 / November 1993
- Published online by Cambridge University Press:
- 26 April 2006, pp. 343-401
-
- Article
- Export citation
-
In this work we consider how nonlinear hydrodynamic effects can lead to a mean force of interaction between two spheres of equal radius a undergoing translational fluctuations parallel or perpendicular to their line of centres. Motivated by amplitudes and Reynolds numbers characteristic of Brownian motion in colloidal systems, nonlinearities due to motion of the boundaries and to inertia throughout the fluid are treated as regular perturbations of the time-dependent Stokes equations. This formulation ultimately leads to a prescription for computing, at leading order, the time-average nonlinear force for the case of pure oscillatory modes – which represents the Fourier decomposition of more general motions. The associated hydrodynamic problems are solved numerically using a least-squares boundary singularity method. Frequency-dependent results over the whole spectrum are presented for a sphere-sphere gap equal to one radius; illustrative calculations are also carried out at other separations. Subsequently we extend the analysis of nonlinear drift to a Langevin equation formulation of the more complex problem of stochastic motion due to thermal fluctuations in the suspending fluid, i.e. Brownian motion. By integrating (numerically) over the spectrum of frequencies, we quantify how the mutual interactions of all translational disturbance modes give rise, on ensemble average, to a stochastic nonlinear force of interaction between the particles. It is particularly interesting that this net interaction – arising from a zero-mean random force – is of O(1) on the Brownian scale kT/a, even though it represents a small O(Re) correction at each frequency of pure oscillations. Finally, we discuss how the presence of stochastic nonlinear drift would lead to non-uniform equilibrium distributions of dilute colloidal suspensions, unless one adds to the random force in the Langevin equation a cancelling non-zero mean component.
Coalescence, torus formation and breakup of sedimenting drops: experiments and computer simulations
- GUNTHER MACHU, WALTER MEILE, LUDWIG C. NITSCHE, UWE SCHAFLINGER
-
- Journal:
- Journal of Fluid Mechanics / Volume 447 / 25 November 2001
- Published online by Cambridge University Press:
- 30 October 2001, pp. 299-336
-
- Article
- Export citation
-
The motion and shape evolution of viscous drops made from a dilute suspension of tiny, spherical glass beads sedimenting in an otherwise quiescent liquid is investigated both experimentally and theoretically for conditions of low Reynolds number. In the (presumed) absence of any significant interfacial tension, the Bond number [Bscr ] = (Δρ)gR2/σ is effectively infinite. The key stages of deformation of single drops and pairs of interacting drops are identified. Of particular interest are (i) the coalescence of two trailing drops, (ii) the subsequent formation of a torus, and (iii) the breakup of the torus into two or more droplets in a repeating cascade. To overcome limitations of the boundary-integral method in tracking highly deformed interfaces and coalescing and dividing drops, we develop a formal analogy between drops of homogeneous liquid and a dilute, uniformly distributed swarm of sedimenting particles, for which only the 1/r far-field hydrodynamic interactions are important. Simple, robust numerical simulations using only swarms of Stokeslets reproduce the main phenomena observed in the classical experiments and in our flow-visualization studies. Detailed particle image velocimetry (PIV) for axisymmetric configurations enable a mechanistic analysis and confirm the theoretical results. We expose the crucial importance of the initial condition – why a single spherical drop does not deform substantially, but a pair of spherical drops, or a bell-shaped drop similar to what is actually formed in the laboratory, does undergo the torus/breakup transformation. The extreme sensitivity of the streamlines to the shape of the ring-like swarm explains why the ring that initially forms in the experiments does not behave like the slender open torus analysed asymptotically by Kojima, Hinch & Acrivos (1984). Essentially all of the phenomena described above can be explained within the realm of Stokes flow, without resort to interfacial tension or inertial effects.
Shear-induced lateral migration of Brownian rigid rods in parabolic channel flow
- Ludwig C. Nitsche, E. J. Hinch
-
- Journal:
- Journal of Fluid Mechanics / Volume 332 / February 1997
- Published online by Cambridge University Press:
- 10 February 1997, pp. 1-21
-
- Article
- Export citation
-
This paper addresses the cross-stream migration of rigid rods undergoing diffusion and advection in parabolic flow between flat plates – a simple model of a polymer that possesses internal (rotational) degrees of freedom for which the probability distribution depends upon the local shear rate. Unequivocal results on the observable concentration profiles across the channel are obtained from a finite–difference solution of the full Fokker–Planck equation in the space of lateral position y and azimuthal angle φ, the polar angle θ being constrained to π/2 for simplicity. Steric confinement and hydrodynamic wall effects, operative within thin boundary layers, are neglected. These calculations indicate that rods should migrate toward the walls. For widely separated rotational and translational timescales asymptotic analysis gives effective transport coefficients for this migration. Based upon angular distributions at arbitrary rotational Péclet number – obtained here by a least–squares collocation method using trigonometric basis functions – accumulation at the walls is confirmed quantitatively by the effective transport coefficients. The results are extended to free rotation using spherical harmonics as the basis functions in the (φ, θ) orientation space. Finally, a critique is given of the traditional thermodynamic arguments for polymer migration as they would apply to purely rotational internal degrees of freedom.